Kinetic study for hopping conduction through DNA molecules 
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Recent experiments indicated that disorder effect in DNA may lead to a transition of the charge 
transport mechanism from band resonant tunnelling to thermal activated hopping. In this letter, 
based on Mott's variable-range hopping theory we present a kinetic study for the charge transport 
properties of DNA molecules. Beyond the conventional argument in large-scale systems, our 
numerical study for finite-size DNA molecules reveals a number of unique features for (i) the I-V 
characteristics, (ii) the temperature and length dependence, and (iii) the transition from conducting 
to insulating behaviors. 
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The fundamental electronic process in DNA molecules 
(DNA electronics) has received great interest in recent 
years. In addition to a large number of indirect opti- 
cal measurements, recent direct electrical measurements 
for charge transport through DNA molecules revealed 
amazingly conflicting transport behaviors, ranging from 
insulatorQ, 0, 13, 3 , Ohmic conductor [ElEHi semicon- 
ductor 8] , to even superconductor % As a consequence, 
the intrinsic charge migration mechanism remains highly 
controversial. 

While quantum transport (in terms of band resonant 
tunnelling) picture 10, 11, 12] was attributed to the ob- 
served conductor/semiconductor behaviors, Mott's vari- 
able range hopping (VRH) theory was proposed to un- 
derstand the temperature dependence of the optically 
measured conductivity of A-DNA |0, It was also 

suggested that it is the disorder effect leading to the in- 
sulating behavior observed in other experiments |^ |3 Q . 
In particular, the experiment by Yoo et al 15] provided 
clear evidence for a polaron hopping mechanism which is 
responsible to the electrical conduction through a DNA 
molecule with length about 20nm and containing iden- 
tical base pairs. In this Letter, we employ an extended 
version of the VRH model to study the hopping conduc- 
tion through DNA molecules. Differing from the con- 
ventional treatment in bulk or large scale systems based 
on a qualitative argument, we base our analysis on di- 
rect numerical simulation for finite systems, which is of 
particular interest in light of recent experiments 0, Il5l ] . 

DNA molecule with random base pair sequence such 
as the A-DNA 0, 0, ^], or, with identical base pairs 
but influenced by complex environment 0, h5| , can be 
treated as a one-dimensional disordered system, where 
the dominant channel for charge migration is a series 
of hops between the localized states. The thermal ac- 
tivated hopping rate between two localized states, say, 
the mth and nth states separated by a distance Rmn, 
can be described as O, k^n = koe-'^"^--"-^^'""/'"'^ , 
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where /cq is the attempt-to-escape rate, the localiza- 
tion length, and Wmn the energy difference of the two 
states. We identify Wmn — ^a/Rmm where A denotes 
the total energy disorder strength of the system, and a is 
the distance between two adjacent base pairs. Two ad- 
ditional remarks in relation to the hopping model to be 
adopted are as follows: (i) Instead of considering only the 
most probable hops as in the standard VRH theory, the 
present work will take into account all possible hops with 
the probabilities described by k,nn- (ii) The localization 
length a is to be influenced by the structural fluctua- 
tions (i.e. the dynamic disorder effect); thus it depends 
on temperature. Following Yu and Song [l^ . we model 
this effect by a = ao + Q;i tanh(r/Td)^, where ao de- 
scribes the static disorder, and the second term is from 
the dynamic structural fluctuations. 

For the electrical transport measurement, at zero bias 
voltage all the localized states are occupied, resulting 
from the hybridization of the individual HOMO states 
of all the base pairs. Switching on the bias voltage, a 
non-equilibrium state is developed, which is described 
kinetically by the rate equation 

jn — (1 fn ) ^ ^ krjim fm ^ ^ (1 fm ) ^mn fn 

m m 

+ C(l-/«)-^r/n. (1) 

Here /„ is the probability of hole occupation. In this 
equation, two types of hopping rates are involved, i.e., 
hopping between localized states in the DNA molecule, 
and hopping between the (localized) molecular states and 
the electrodes. The former has been given by the stan- 
dard VRH model. In the following we develop an expres- 
sion for the latter, which is characterized by the rates k]^ 
and fc™*. 

In contrast to quantum transport, the classical hop- 
ping considered here involves real transition between the 
electrode and the localized molecular states with differ- 
ent energies, and the individual excess energy is gained 
from or lost into the surrounding environment. In gen- 
eral, this inelastic hopping process is described by the 



2 



non-radiative transition, with rate 
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i7„p|2(4^AfcBT)-5exp 
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Here, 



X is the reorganization energy of the environ- 
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ment, and Epn = Cp — Cn, denoting the energy dif- 
ference between the electrode state with momentum p 
and the nth locahzed molecular state, which are coupled 
with strength iJ„p. Physically, Hnp — H exp{—f3Rn/2), 
where is the tunnelling length which is assumed 
here independent of the electrode states, and i?„ is the 
distance between the electrode and the nth localized 
state. Also, in later numerical evaluation of the hopping 
rates between the molecular states and the electrodes, we 
would approximate the energy difference Epn ^ ep — eo, 
where eo denotes the energy of the localized molecular 
state nearby the electrode. This approximation makes 
sense in viewing the rapid decay of 7Jp„ with The 
total hopping rate from the electrode to the localized 
molecular state can be evaluated by integrating the elec- 
trode states. For instance, the hopping rate from the 
left electrode to the nth molecule state reads, fc„i, = 
gL J dtphnpfhi^p — Ml), where gL is the density of states 
of the (left) electrode. In the following numerical cal- 
culation, the combined parameters ^l{R) — '^T^\H\'^gL{R) 
will be used to characterize the coupling strength be- 
tween the molecule and electrode, and will be commonly 
adopted as 0.2 meV. The Fermi level of the left electrode 
in the presence of applied voltage V is assumed to be 
^J'L = Ego + eV/2, where the gap energy Ego = Ep - eo, 
with Ep the equilibrium Fermi energy of the electrode. 
Other hopping rates between the electrode and the local- 
ized molecular states, i.e., k^n, f^nn and knB.i can be sim- 
ilarly evaluated. Thus the rates in Eq. |^ are obtained 

as fc™ = knL + knR, and = ^Ln + kpjn- 

After identifying all the rates in Eq. J^l, we can obtain 
the time-dependent evolution of the occupation proba- 
bilities on the individual localized states in response to 
an applied voltage. In this work we are in particular in- 
terested in the stationary hopping current through the 
DNA molecule, which can be evaluated as 



I In)knL - kLnfn], 
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under the condition /„ = 0. 

Figure 1 shows the I-V characteristics associated with 
the hopping conduction, where the DNA molecule with 
30 base pairs is exemplified. Here the calculated cur- 
rent is eventually saturated at certain bias voltage, as 
shown in Fig. 1(a), owing to adoption of the simplified 
one-band model. In this work we will use the saturated 
maximum current (/max) to characterize the conduction 
ability (equivalent to the average conductivity over dif- 
ferent voltages). The voltage gap in Fig. 1(a) is roughly 
determined by the relative position of the HOMO level of 
the DNA base pair near the electrode from the Fermi sur- 
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FIG. 1: I-V characteristics under hopping conduction 
through DNA molecule with A'^ = 30 base pairs. Plotted are 
results for, (a) different energy gaps at temperature T = 300 
K, and (b) different temperatures with a given Ego — 0.3 eV. 
Other parameters adopted here are the reorganization energy 
A = 0.1 eV, and the disorder energy A = 0.15 eV. 



Note that this gap differs from its counterpart in the 
quantum transport regime, where the individual base- 
pair states interact with each other and an energy band is 
formed, and the voltage gap is determined by the distance 
of the upper e dge of the energy band from the electrode 
Fermi surface \ldL IllL IT^. In different experiments, this 
gap may differ considerably, leading to cither the metallic 
ohmic or the semiconductor behaviors. We thus adopted 
several values of Ego in Fig. 1(a) to illustrate the possi- 
bly observed I-V characteristics. Moreover, the hopping 
conduction displays a characteristic temperature depen- 
dence of thermal activation, as shown in Fig. 1(b), which 
is in good agreement with the experiment by Yoo et al 
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FIG. 2: Maximum (saturated) current versus inverse temper- 
ature. The adopted parameters are A = 0.1 eV, Ego = 0.3 eV 
and = 30. The inset shows the thermal activation energy 
Ea versus the disorder strength A. 

The thermal activation characteristics are further man- 
ifested clearly by the exact exponential dependence of the 
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inverse temperature at high temperature shown in Fig. 
2, where the slope gives the thermal activation energy 
Ea- In general, the thermal activation energy depends 
on the disorder, as quantitatively shown in the inset of 
Fig. 2. As an illustration, for disorder A — 0.2eV, we 
obtain Ea = 0.12eV, which agrees well with both the 
experiment [T^j l and the ab initio calculation W^. Low- 
ering the temperature, there appears a notable regime in 
which the conduction is of weak dependence of the tem- 
peratures. The transition takes place at temperature of 
200 ^ 250K, which is again in agreement with the exper- 
iment |l5j| . The results numerically obtained here can be 
qualitatively understood by the VRH argument 0|. 
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FIG. 3; Ohmic behavior of length dependence under hopping 
conduction. The plotted lines correspond to different disorder 
strengths. Other parameters are A = 0.1 eV, Ego = 0.3 eV, 
and T — 300 K. The inset displays the effective resistivity 
(obtained from the slope of the lines in the main part of the 
figure) versus the disorder strength, and the conducting-to- 
insulating transition (i.e. with current from nA to pA) by 
increasing the disorder for a DNA with A = 30 base pairs. 

Another characteristic feature associated with the hop- 



ping conduction is the ohmic behavior of length depen- 
dence, as shown in Fig. 3. This feature differs from either 
the coherent tunnelling through a disordered system at 
zero temperature, or the quantum transport through sys- 
tem without disorder: the former has the characteristic 
length dependence (e^^"^ — 1) while the latter 
leads to a maximum current contributed from the en- 
tire HOMO band which is almost independent of the 
molecule length. The disorder strength would signifi- 
cantly affect the conduction property as manifested in 
the inset of Fig. 3, by the effective resistivity (i.e. the 
slope) and the conducting-to-insulating transition (i.e. 
with current from nA to pA) by increasing the disorder. 
The insulating transition also happens by increasing the 
molecule length. As a rough estimate, consider the hop- 
ping conduction through DNA molecule at room temper- 
ature and with energy disorder A — 0.15 eV. From the 
result in Fig. 3, we obtain an estimate for the maximum 
current which would decrease from 60 pA to 0.5 pA as 
the base-pair numbers increase from 30 to 3000, i.e., to 
the length of micron. This is nothing but the insulating 
transition of DNA molecules on micron scales, which has 
been commonly concluded in a number of recent experi- 
ments 0,1111. 

In summary, we have presented a kinetic study for the 
transport properties of disordered DNA molecules, based 
on Mott's variable-range hopping theory. A number of 
unique features associated with the thermal activated 
hopping mechanism were discussed with respect to either 
the already known experimental results, or the possible 
future experiments. 
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